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Abstract 

In this paper, we present two adaptive methods for the basis enrichment of the mixed Generalized 
Multiscale Finite Element Method (GMsFEM) for solving the flow problem in heterogeneous media. We 
develop an a-posteriori error indicator which depends on the norm of a local residual operator. Based on 
this indicator, we construct an offline adaptive method to increase the number of basis functions locally 
in coarse regions with large local residuals. We also develop an online adaptive method which iteratively 
enriches the function space by adding new functions computed based on the residual of the previous 
solution and special minimum energy snapshots. We show theoretically and numerically the convergence 
of the two methods. The online method is, in general, better than the offline method as the online 
method is able to capture distant effects (at a cost of online computations), and both methods have 
faster convergence than a uniform enrichment. Analysis shows that the online method should start with 
certain number of initial basis functions in order to have the best performance. The numerical results 
conhrm this and show further that with correct selection of initial basis functions, the convergence of 
the online method can be independent of the contrast of the medium. We consider cases with both very 
high and very low conducting inclusions and channels in our numerical experiments. 


1 Introduction 

Many real-world problems involve multiple scales and high contrast. To solve these problems, we often adopt 
some forms of model reduction such as upscaling and multiscale methods. These methods can reduce the 
degrees of freedom of a problem. For example, in upscaling methods [M113 mi HU, the multiscale media 
are upscaled so that the problem can be solved on a coarse grid. In multiscale methods [S1E112D11I711IHIII1 
ESlHIHliniEllIllEaiBEH], basis functions are solved on a fine grid to capture the multiscale features of a 
medium and the problem is then solved on the coarse grid with these basis functions. 

In this paper, we will present two adaptive enrichment algorithms for the generalized multiscale finite 
element method (GMsFEM) in solving the mixed framework of the flow problem in heterogeneous media 
[^. The first method is based on a local error indicator. We use this indicator to search for the regions, 
where more basis functions are needed. This method will only add pre-computed basis functions, which are 
computed in the offline stage so we call it an offline adaptive method. In the second method, new basis 
functions are computed based on the previous solutions. We call it an online adaptive method. 

GMsFEM is a generalization of the classical multiscale finite element method [35]. In the classical method, 
one basis function per coarse edge is used to capture the multiscale features. For the multiscale mixed hnite 
element method, one may see |1I3[I|. GMsFEM allows more basis functions per coarse edge to be used to 
take into account the effects of non-separable scales. The main idea is to solve local spectral problems for 
the selection of basis functions. The formation of basis functions in GMsEEM can be divided into offline 
and online stages. In the offline stage, offline basis functions are computed based on the multiscale features 
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so that these functions can be reused for any input parameters to solve the equation. Online functions 
are those depending on the parameters. In an adaptive algorithm is developed to enrich the space by 
adding basis functions which are formed in the offline stage. In m, adaptive methods which involve the 
formation of new online basis functions based on the previous solution are developed. These methods show 
significant acceleration in the convergence rate of GMsFEM. There are also related methods developed for 
the discontinuous Galerkin formulation in [10] and [12] . 

In the paper, we will focus on the mixed framework of the flow problem. The mixed methods are important 
for many applications, such as flows in porous media, where the mass conservation is essential. We developed 
two adaptive methods to enrich the function space. One involves only offline basis functions while the other 
adds new online basis functions that are constructing using special minimum energy snapshots. We call 
them an offline and an online adaptive methods respectively. Two local spectral problems are developed for 
constructing multiscale basis functions. Both of them can be used in the online method, but only one can be 
used in the offline method. We propose error indicators which are based on the and the H{div) norms of 
the local residual. These error indicators can be used to approximate the error of the solution. From [5], we 
know that the error between the GMsFEM solution and the fine grid solution involves two parts: one due 
to the selection of the basis functions and the other due to the discretization of the source function. In this 
paper, we will assume the error due to the discretization of the source function is small and consider only 
the former part. The offline adaptive method depends on the error indicator to help the selection of basis 
functions. The online adaptive method produces new basis functions iteratively by projecting the previous 
solution on the space of divergence free functions. We emhasize that [5] gives a-priori error estimate of the 
mixed GMsFEM, and the purposes of this paper are a-posteriori error estimates and adaptivity. 

In our analysis, we prove the convergence of the two methods. It can be shown that the error is bounded 
by the local error indicators. By adding offline basis functions to those coarse grid edges with large error 
indicator, we can show a guaranteed convergence rate for the error of the solution together with the local 
error indicators. The convergence rate depends on the parameters of the offline adaptive method. These 
parameters control the number of coarse grid edges to be chosen and the number of basis functions to be 
added for each of those edges. For the online adaptive method, a set of non-overlapping subsets of the domain 
is selected. New basis functions are computed on each of these subsets. We show that the convergence rate 
depends on the norm of the residual operator restricted on those subsets, and also the eigenvalues of the 
offline basis functions that are not included in the initial basis. 

We present some numerical results to show the convergence behaviour and some properties of the adaptive 
methods. We consider both high and low conductivity inclusions and channels in the domain. By comparing 
the adaptive methods to uniformly enriching the function space, one can see the efficiency of the adaptive 
methods. In particular, the performance of the online adaptive method is generally the best since it adds 
functions which are computed based on the previous solution while the offline adaptive method enriches 
the space by adding basis functions which are independent of the input parameters. We will see that both 
the choice of the non-overlapping regions and the initial number of basis functions on each coarse grid edge 
affect the convergence of the online adaptive method. Some of the eigenvalues from the spectral problems 
are sensitive to the contrast of the problem. By including functions corresponding to those eigenvalues in 
the initial basis, the convergence of online adaptive method becomes independent of the contrast. 

The rest of the paper is organized in the following way. In the next section, we briefly introduce the 
basic idea of mixed GMsFEM. At the end of the section, we give the detail of the two adaptive methods. 
In Sectionwe state and prove the convergence results for the adaptive methods. In Section]^ numerical 
results are given to illustrate the convegence behaviour of the adaptive method and the factors affecting the 
convergence. The paper ends with a conclusion. 
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2 Method description 

2.1 Overview 

Consider the high-contrast flow problem in a mixed formulation: 

K~^v + \/p = 0 in D 
div(v) = / in D, 

with Neumann boundary condition v • n = g on dD, where «: is a high-contrast permeability field, D is the 
computational domain in M" and n is the unit outward normal vector of the boundary of D. 

We will solve the equation on two meshes with different scales. Let be a partition of D into finite 
elements (triangles, quadrilaterals, tetrahedra, etc.), where H is the mesh size. We call the coarse grid. 
Next we construct a finer grid. For each coarse grid element K G , we further partition K into a finer 
mesh such that the resulting partition of D with size h is conforming across coarse-grid edges. We call 
the fine grid. Denote the set of all faces of the coarse grid as 5^, and let be the total number of faces 
of the coarse grid. We define the coarse grid neighborhood uji of a face Ei G as 

LOi = \J{K gT^ :E,G dK}, 

which is indeed a union of two coarse grid blocks. 

Next, we define the notations for the solution spaces for pressure and velocity. Let Q be the space of 
functions which are constant on each coarse grid block. We will use this space to approximate p. For the 
velocity space, we will first construct a set of basis functions /3inap for each coarse grid neighborhood 
We call 14nap = ©BiGfTH L^nap the snapshot space, where Fsnap = span ^/Sinap^. The snapshot space is an 
extensive set of functions which can be used to approximate the solution v. However, this space is large and 
we will reduce it to a smaller one before we solve the equation. From each t4nap, we select a set of basis 
functions /3mL Denote Hmi = span ^/3ms^ and 14is = V^s ■ The size of Fms is generally smaller than 

14nap- We will use the space Hms to approximate the velocity v. 

With the pressure space Q and the velocity space Fms, we solve for € Q and Vms € Fms such that 



/ div(w)pms = 0 
Jd 


yw G 


with boundary condition v^s • n 
projection of g in the sense that 



= [ fq Vq e Q, 

JD 

= {u e Fins : V ■ n = 


( 2 ) 


0 on dD}, and gn is the 


/ {gn - g)(l) ■ n = 0 y(j)G ^^np and E^ C dD, 

JEi 

and gn is constant on each fine grid face. 

For n C D and v G Fsnap, we define two norms ||u||c 2 (a;K-i) and ||u||B(div;a;K-i) as 

lkll£2(n;K-i) = 

and ^ 

lkllff(div;n;K-i) = f [ div{v)A . 

\Jn Jn J 

We will use these two norms in the adaptive methods. 

In the coming sections, we will discuss the formation of the snapshot space I4nap and the method of 
selecting the /3m^s. We will also give two adaptive methods of enrichment of the multiscale space Fms so as 
to get a more accurate solution without using too many basis functions. 
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2.2 Snapshot space 

In this section, we will present the construction of the snapshot space which is a large function space 
containing basis functions up to the resolution of fine grid faces on the coarse grid faces. We construct 

the local snapshot bases /3snap by solving a local problem on each coarse grid neighbourhood uii, and then 
(i) 

combine the ^snap’s to form the snapshot space Vijnap- 

Let Ei G 8^. We write Ei = u;=i Cj, where the e^ ’s are the fine grid faces contained in Ei and Ji is the 
total number of those fine grid faces. We will solve the following local problem to obtain /3snapj 

n~^vf + \/pf = 0 

div(Uj-*^) = 

subject to the homogeneous Neumann boundary condition = 0 on duJi. We want the snapshot basis 

to contain solutions of the local problem with all possible boundary conditions on the edge Ei up to the fine 
grid resolution. Therefore, the problem is solved separately on each coarse-grid element K C uji with the 
additional boundary condition up) • rii = on Ei, where <5^*^ is defined by 

^ [0 on Ei\ej, 

(0 

and rii is a fixed unit-normal vector for each face Ei. The function ' is constant on each coarse grid block 
and it should satisfy the condition ^ • rii for every K C uii. 

The set of solutions to ^ is the local snapshot basis ,dinap. Using [Je-^£h Psnap as a basis, we have the 
snapshot space Usnap- 


in uji, 
in uji, 


(3) 


2.3 Offline space 

After we obtain the snapshot spaces I4nap, we perform a dimension reduction to get a smaller space. Such 

reduced space is called the offline space. The reduction is achieved by solving a local spectral problem on 

(i) 

each coarse grid neighborhood Wi. Some of the eigenfunctions will be used to form the local basis /3ms. The 
local spectral problem is to find real number A and v G V such that 

a{v,w) = Xs{v,w), yw G V, (4) 

where V is the snapshot space, a and s are symmetric positive definite bilinear operators defined on U x U. 

We propose the following two possible spectral problems for the basis selection. In these two problems, 
V is set to be I4nap. 


Spectral problem 1: We take 


ai{v,w) 


K ^{v ■ ni){w Hi), 


Si{v, w) 



div(u)div(r(;) 


for V, w G I4nap, where rii is a fixed unit normal vector on Ei. The term 1/iL is added so that Ui and Si 
have the same scale. 
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Spectral problem 2: For v € Fsnip, we define v to be the extension oi v ■ rii in V's^nip by 

minimizing the energy norm on uji, i.e., we find v € ^"^^p such that v ■ rii = v ■ rii on Ei and 

l|w||£ 2 (^.^-i) < for all ip e We take 


a^(y^yj)= f K-^V-W, 
J UJi 

S^{v^w)= j K-^-W. 

J UJi 


Remark that the eigenvalues of this spectral problem is always bounded above by 1. 


(■i) 

Note that spectral problem 1 involves only the space Fsnap- However, in spectral problem 2, the calculation 
of V involves all the Fsnap’s such that uij and uJi have non-empty intersection, and so spectral problem 2 
requires slightly more calculation. We note that v is the minimum energy extension oi v ■ ni\Ei in the space 

®ujnLJiP<p f^snap. 

After solving the spectral problem in a coarse grid neighborhood Wj, we arrange the eigenfunctions 
in ascending order of the corresponding eigenvalues 

(i) 

We then let /3ms be the set of the first li eigenfunctions, where U is a predefined integer or li is depending 
on the eigenvalues. The selection of li will be discussed in the analysis section. When the spectral problem 
is specified and the k's are selected, we construct the offline space Vms- 

Suppose li is fixed for each coarse grid neighborhood oJi. We define Amin = min^.g^ff The value of 

(Amin)”^ indicates the error between the multiscale solution and the solution obtained by using the whole 
snapshot space. Therefore, we want (Amin)”^ to be as small as possible. However, in spectral problem 2, 
we can see that (Amin)”^ is bounded below by 1. Therefore, this spectral problem will only be used in the 
online adaptive method, and is shown to be crucial in selecting initial bases. 


2.4 Offline adaptive method 

In this section, we will introduce an error indicator on each coarse grid neighborhood. Based on this estimator, 
we develop an offline adaptive enrichment method to solve equation Q iteratively by adding offline basis 
functions supported on some coarse grid neighborhoods in each iteration. In this offline adaptive method, 
we will use spectral problem 1. 

For each coarse grid neighborhood Wi, we define the residual operator Ri as a linear functional on Fsnap 

by 

Ri{v) = / K~'^Vms-V- 

J UJi 

where (umsj Pms) is the solution pair of ([^. 

We take as our error indicator, where 

1. ^snap ) 


div(u)pn 




y{i) 

* snap’ 


\R,\ 


K, 




= sup ■ 
4nlp 




«6v;V 


\H(di 


IV.iOi'.K 


This quantity can be used to approximate the error since we have 

ll'^^snap ~ < Cerr ^ ^ ll^illyU) + ; 
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where fsnap € V"snap IS the solution obtained by solving the equation using the whole snapshot space 14 nap- 
We will prove this inequality in the next section. 

With this error indicator, we now present the offline adaptive method. In this method, only spectral 

problem 1 will be used. We let m > 0 denote the enrichment level. For each coarse grid neighborhood Wi, 

(i) 

we use to represent the number of eigenfunctions used to form V,is ■ We use to denote the space Fms 
at enrichment level m. 


Offline adaptive method: Fix the number 6 and Sq with 0 < 6 *, Jg < 1. We start with iteration number 
m = 0. Fix an initial number of offline basis functions for each coarse grid neighborhood to form the 
offline space Then, we go to step 1 below. 

Step 1: Find the multiscale solution. We solve for the multiscale solution and G Q satisfying 


K 


ID 


[ divHp-=0 

JD 

[ div«g)g= f fq 

JD JD 


yw G (F™)° 


yq G Q. 


(5) 


Step 2: Compute the error indicators. For each coarse grid neighborhood Wj, we compute the local residual 
r]f = and rearrange the residuals in decreasing order r/i > 772 > • • • > rjN^- 

Step 3: Select the coarse grid neighborhoods where basis enrichment is needed. We take the smallest k 
such that 

Afe fc 

i=l 


We will enrich the offline space by adding basis functions which are supported in the coarse grid 
neighborhoods wi,..., 


Step 4: Add basis functions to the space. For each of the selected coarse grid neighborhood uii from step 
3, we will take the smallest Si such that We then set = I™ + Si so that Si 

more basis functions are included to form a larger space For the other neighborhoods, we set 

— ^T- We will explain the reason for choosing such Si in the analysis section. 


After Step 4, we repeat from Step 1 until the global error indicator 'Hi is small enough or the 

total number of bases reaches certain level. The calculations of all the local error indicators can be time 
consuming. However, since the error indicators are independent of each other, the computation can be done 
in a parallel approach. 


2.5 Online adaptive method 

Next, we will present another enrichment algorithm which requires the formation of new basis functions 
based on the solution of the previous enrichment level. We call these functions online basis functions as 
these basis functions are computed in the online stage of computations. With the addition of the online 
basis functions, we can get a much faster convergence rate than the offline adaptive method. 

We first define a linear functional which generalizes the residual operator in the offline adaptive method. 
Given a region fl C D, let Vq be the space of functions in Fsnap which are supported in H, i.e. Vq = 
0^.ca A^snap- Let Vii denote the divergence free subspace of Fq. Define the linear functional Rq, on Fn by 

Rn{v) = [ -v - [ div(u)p™g Vu G Fq. 

JD JD 
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The norm ||i?a||y^, we use, is given by either 


\\Rn\\v^ = sup ' 

vGVn 


Rniv) 


kllff(db 


or \\Rn\\v^ = sup ■ 
vGVa 


Rail 


kll£2(n;K-i 


depending on which spectral problem is used. Remark that if we restrict Rq on Vh, we have 


Rn{v) = [ K -v Vu e Vq. 

JQ 

For the case of id = Wi for some i, Rn is the same as the residual operator Ri in the offline adaptive method. 

Similar to the offline adaptive method, we use m to indicate the enrichment level and to denote the 
velocity space at enrichment level m. In the online adaptive method, we can use either spectral problem 1 
or 2. However, since the online basis functions constructed in each enrichment level are divergence free, we 
must ensure that before any enrichment of the space, there is at least one basis function supported on each 
coarse grid neighborhood, which is not divergence free. 


Online adaptive method: Let m = 0. We start by choosing an initial number of offline basis functions, k, 
for each coarse grid neighborhood We use the first li eigenfunctions from each coarse grid neighborhood 
uJi to form the initial velocity space If the first li eigenfunctions of a coarse grid neighborhood are all 
divergence free, we may artifically construct an extra basis function that is not divergence free and include 
it into We go to step I below. 

Step 1: Find the multiscale solution. We solve for the multiscale and G Q as in step I of the 

offline adaptive method. 

Step 2 : Select non-overlapping regions. We pick non-overlapping regions fli, il 2 , •. •, C D such that 
each U,j is a union of some coarse grid neighborhoods. 

Step 3: Solve for online basis functions. For each flj, we solve for (pj £ Vq,. such that 

Rq^{v) = -v Vn £ Vq., 

jQ.j 

i.e., we solve for the Riesz representation of Rq- in Vfi .. Those (j)j^s are the new online basis functions. 
We update the velocity space by setting F,^+^ = 0 span{(/)i, (j) 2 , •. ■, (pj}- 

Again, after Step 3, we repeat from Step 1 until the global error indicator is small or we have certain 
number of basis functions. 

In our calculation, after obtaining the online basis function pj in step 3 of the method, we sometimes 
normalize it before computing the matrix in the finite element method. We can see that pj is a projection 
of the multiscale solution on the space Vq.. We have \\Pj\\c'^{D-K-^) = ■ When is close 

to the snapshot solution (the solution Vms of equation (0 with Fms being the whole snapshot space) in the 
region by the first equation of ([^, the norm of the projection pj will be small. Adding pj directly into 
calculation will make the matrix in the calculation close to singular. 


3 Convergence analysis 

In this section, we will present the proofs for the convergence of both the offline and the online adaptive 
method. We first define some notations that will appear in the results. 
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We denote the maximum number of faces of a coarse grid block by Nq-. Let Usnap G 14nap and Psnap G Q 
denote the snapshot solution, i.e. 'Cgnap and Psnap satisfy 


ID 


K ^Vsni^p ■ W 


- / div(w)psnap = 0 Vw G V;°ap, 
J D 

/ div('Usnap)g = / fq ^q&Q, 

Jd Jd 


( 6 ) 


with Usnap ■ n = gu on dD. 

With these notations, we have the following result for the error indicator. 


Lemma 1. We have 


We 


W'^snap '^ms 




If spectral problem 1 is used, then Cerr = and || • || (o is || • \\Hldivu> -K,-^)- If spectral problem 2 is used, 
then Cerr = Nj- and |j • ||„(i) is || • The value of Cv depends on the polynomial order of the fine 

grid basis functions in Vsnap ■ 

Proof. For any v G 14napj we can see from the construction of V"snap that div(?;) is constant on each coarse 
grid block. Hence, by the second equation of © and § , we have 


I div('Usnap "Cnis) — / div(Usnap)div('Usnap "Cms) 

Jd Jd 

- / div(?;ms)div(Wsnap - Vms) 

Jd 


= 0 . 


Thus, div(Ums - '^snap) = 0. 

Next, since Ugnap - ^'ms G V^nap, we have 


^ I'Csnap 'amsi — 


f K-\ven.p-V^e)-{Vsn.p-Vms) 

Jd 

/ div(Usnap ^ms)(Psnap Pms) 

Jd 

Using the first equation of Q, we get 

/ «-'ks„ap-^msP=- / K-\^e-iVsn.p-V^e)+ f div(^;anap - 

Jd Jd Jd 


By definition, we can write 


I 1^ l^snap ^ms| — {Rd , ^snap ^ms)- 

Jd 


We decompose r^snap ~ as the sum of functions from the Usnap’s, i.e. t'snap ~ Vms = where 

fh) G i4nap. Each rih) can further be written as sum of two components: one in Umi and the other in 
span(/3inap \ /?ms). Let Vt''^ be the latter one. We get 

We 

^Rd-, 'Cgnap ’^^ms) — ^ '. {Ri^ V^ '^) 
i=l 
We 



Using the definition of the spectral problems, we get 


N, 


N, 


/(O 


2 = 1 


2=1 

N, 


2=1 


where || • Hv'^. = || • ||//(div;a;i;K-i) if spectral problem 1 is used, and | if spectral problem 

2 is used. 

Next, we consider the two spectral problems separately. 


Spectral problem 1: For each i, we have 

[ K~^ {{Vsnap - Vms) ■ riif. 

' JEi 

Thus, by Cauchy Schwarz inequality. 


{Rdi 




2=1 
iVe 




< 


< 




'Ei 
1 

2 / iVe 


\2=1 


CvH 


mlyi. y 


' N, 


^ ((“^snap ^ms) ‘ ^2) 


^2 = 1 


'y ' II nAi) j Il'T’snap 


\i=l 


where Cy is a constant depending on the polynomial order of the fine grid basis functions in 14nap- 
Spectral problem 2: Similar to spectral problem 1, for each j, we have 

J (jJi 

< / 

J UJi 

where we used the minimum energy property of 7 ;b). Therefore, 


^ I'^^snaD 


snap *^ms | 5 


{Rd. ^^snap - ^^ms) < ^ H (V «p)* )” M (E ^ 


^ l^snap ^ms| 


' iVe 


< \JN'Y 11^^’bdiL)* J ll'*^snap ■L'ms||£2(^.,..-l). 


□ 


From the proof of Lemma the error between the multiscale solution and the snapshot solution is 
bounded above by the norm of the global residual operator Rd, which in turn can be estimated by the sum 
of the error indicator ||ii.i||, (i) ...(A^y )“^. 

I, » snap / ^ 
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Next, before we show the convergence of the offline adaptive method, we need a local version of the 
inf-sup condition in [^. We will show the proof of this simplier case and compute the constant in the result. 
Here is the statement. 


Lemma 2. For coarse grid neighborhood uji, write LOi = Ki U K 2 where Ki and K 2 are the two coarse grid 
blocks composing uJi. Then, for any p G Q, we have 


\\p - pWc^oj.) = sup 


J div{v)p 


where p = p and Cl^^p is the infimum of over all v G Vml with v -Ui = 1. 

Proof. Note that 

/ div(ti)p = p V ■ n = 0. 

J u)i J duji 

We may assume p = 0. Let po = p/\\p\\c'^(uj^). Using this notation, we have 

/ div(u)p= ||p||£2(^_) (poki / div(u)-hpolifs / div(u) 

Juji \ JKi JK2 

= \\p\\c^iu,,)iPo\Ki-Po\K2) 

JEi 

Finally, we evaluate (polifi ~ Po|if 2 )' Since p^ = 0, we have 

l^ll(Polifi) + |-f^2|(po|if2) = 0- 

Using |bo||£ 2 (,^,) = 1, we get 


\Ki\{p[j\Ki) + \K2\(p[j\K2) — 1- 


Using these two, one can check that 


Hence, we have 


Po\ki — Po\k 2 — 
L, div(u)p 




1 

m 


sup 


= \\p\\cHuJi)] 


Kl\ ' |if2| |b||£2(c^i;K-i) ’ 


Ie, 


which completes the proof. 


□ 


We define some symbols before we move on to the proof of the convergence. Let i?™ denote the residual 
operator Ri using the solution p^s)- We define 


sr = WK 


I 

V * snap 


.(A« .)- 


(7) 


This symbol S'™ is indeed the error indicator pi at the m-th enrichment level. We have the following lemma 
for this symbol. 


Lemma 3. Let be the expression defined in 0. Then, for any a> Q, we have 


{ST+^Y<[l + a) 


'ip+i 


(d 




m+1 


" + 1 


2 


where = 2(a|^+i_|_^) ^(max{C®^, 1})^ and Cl’fff is the constant from Lemmaj^ at the m-th enrichment 
level. 
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Proof. For any v € Vs 


(i) 


snapi 


using the definition of we have 


RT+\v) = RT{v) + [ -vZ)-v- [ div(z;)(p-+i - p™ ). 

J iOi UJi 

Taking supremum with respect to v and noting that T4nap^^ = f4nap, we get 


s: 


m+1 


< 


\(i) 

V+' + l, 


sr+i, 


where 


Next, we have 


/=(A«„^^)t sup 


L, ^ - <s) ■ V - L div(n)(p™+i - p™ ) 


vevjalp 


lH(div;uii;K-^) 


[ K - v^J -v- f div(?;)(p™+i - p™ ) 

J UJi ^ UJi 

-- [ - <s) -v- f div{v){p^+^ - p™ - (pZ+^ - pSJ) 

J IjJ'i 'J UJo 


< - <Jlc2(a;.;«-i)lkll£2(p;.;«-i) + I! div(n) || £2 ||p™+l - p™ - (pZ+^ - P™ ) || £2 

where (pm^^ — pSs) is the average value of — p^s over oji. Using Lemmawe get 


- pZ - - PZ)\\cH.d < c:z sup 


^ ^sup 


< Cl, 


sup 

/^div(i;)(p™+^ -pZ) 

- lbll£2(<^,;K-l) 

vevllJ’' 

sup 

/d - <s) • 

- |bl|£2(p;,;K-l) 

vevllJ’' 

II ‘-^ms 



where Cl;l^ and Umi’™ denote the constant Cs^p in Lemma and the space Umi 
Hence, we estimate I as 


at enrichment level m. 


Therefore we have 


/<(A;ip^^)-5b2max{C:;; 


t ’ -‘-ill ^ms 


£2(0;^;^-!). 


S-^+i < 



.Sr + (A|l.+,)-^v^max{C:, 


sup 






And so, we get 




< (1 + a) 


\ (d 
^ir+i 


A. 


{srr + {i+a-^)Di, 


I m+1 _ m I 
I‘-^ms ‘-^msl 


£2(^-^-l) 5 


^+1 


where Dl^ = 2(A;j^+p^^)-i(max{Q-™, l})^. 

Using this lemma, we have the following result for the convergence of the offline adaptive method. 


□ 
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Theorem 4. Using the notations in the offline adaptive method, there exist positive constants Sq, p and a 
decreasing sequence of positive numbers {Lj} such that the following contracting property holds 


WVsnar, - ( \\Vsnar, - || 

i=l \ 

for any j < m, where Jq o-n-d p satisfy 

<So<l-il- 6o)0^ < P < 1, 

V+i+i 


N, 






2=1 


for any coarse grid neighborhood oJi that are selected to add basis functions, and 

_ CerrLj “t“ p 

CefrLj 1 

The definition of {Lj} is given by (^. 

Proof. In the offline adaptive method, we fixed 0 < ^ < 1 and we choose an index set I such that 








2=1 


iGl 


We write 

+E(^r+')^- 

2=1 

Using Lemma we have 


iVe 

E 

i=l 


\iU 

ip + 1 


<E (i+«)UE^(5rr + (i+«-')i?; 


II m+l _ m I 
II ^ms *^11181 


2 


iGl 


" + 1 


+E ((1+ c^)isrr + (1+ a-^)DUvzp - 


We define 

Lm = NrV + a~^) max 

Ei^Sn 


( 8 ) 


and we assume the number of additional offline basis functions in each enrichment level is chosen such that 


max 

iG/ 





+1 


< ^0 < 1, 


where i5o is a fixed constant. Using this, we get 


E(5r'r <(!+«) E(^rr 


2=1 


2=1 


iVe 

(1 + a)(l - <5o)0^ E(^™)' + 
2 = 1 


We let p = (1 + a)(l — (1 — (5o)6*^) and take a small enough so that 0 < p < 1. Observed that {L^} is a 
decreasing sequence, so we may take any j < m. We now have 

2=1 2=1 
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Note that div(z)™g+^ — = 0 and . Therefore, by the first equation of ([^ and ([^, we 

get 


D 


f^snap 


- - U™ 1 = 0 

‘^ms ) V*^ms ^ms/ 


and so 


which means 


\\v 


snap 


2 

/:2(D;k-1) 


— It'^snap 


— V 


m+1 1|2 

ms \\C^{D\k- 


1) + ||u 


m+1 

ms 


2 


m+1 


2 

£2(£>;«-1) 


||u 


snap 


2 

Z:2(D;«-1) 


||u 


snap 


— V 


m+1 II2 

ms \\C^{D-k-'^)- 


Putting this into (§, we get 


N, 


^^snap 


m+1 II2 


ms WC^iD-K-^) ^ ^ 




N, 

llUsnan " 

i=l 


Csnap — Cms||£2(£,.„-l) -f — 


From Lemma we have 


Afe 


I Wsnap 


\cHD-,.-^)<Ce„J2isrf 


Hence, 


1 

ksnap - + T ^ - 

z=l 


Psnap 'L’ms 


|£2(D;k- 1) + (/3C'err + y) ^(5'™)^- 


We 


i=l 


Finally, we take f3 = to get 


We 


Iksnap - <+' 




'?^snap 1'mslli;2(_D;K-i) 


,m II2 


1 -, 




i=l 






i=l 


which completes the proof. □ 

From Theorem 1^ we can see that the convergence rate depends on the two constants 9 and Jqj which 
are fixed before we carry out the enrichment algorithm. The constant 9 controls the number of coarse 
grid neighborhoods, where enrichment is needed. And the constant Sq is related to the number of basis 
functions we have to add in each coarse grid neighborhood. Note that we have the following inequality for 
the convergence rate 

e >l-ilE£)E. 

Ceri-Lj + 1 

Thus, in order to have a fast convergence, we need a small i5o and a large d, which means there is a tradeoff 
between the convergence rate and the number of basis functions used. 

Now, we state the convergence result for the online adaptive method. 


Theorem 5. Using the notations in the online adaptive methods. We have 


\Vsnap-V^+^\ 


2 

£2(r);, 




I], = l 


1 - . 


C, 




R, 


'v; 


{i) 


(snap’' 




-1 


^snap 


2 

£2(£);K-1)- 
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(10) 


Proof. For any v € with div(t?) = 0, we have 


I V — 

\ t^snap ^ms 


' '^\\c^{D;k-^) — Iksnap 


+ 2 K (Vsnap - ^ 


ID 
f^snap 


— Ill'sna.n Wc^iD^K-^) 

m+l ||2 


|£2(£).k-1) 


\C^{D-k-^) 


> ll'^snap 'I'ms II£ 2 (£).k- 1 ) 


Let the new online basis functions 4>l^ </' 2 , ■■■,4’J be normalized such that \\4'j\\c'^{Q.j-,K-'^) = 1- Let v = 
+ a\(j)i + • • • + ajcjij where aj = K~^v^g ■ (j)j = ||y^ . Check that v is divergence free. 

Using ([To|, 


l^'snap ||£2(£3.„-l) < Ur'snap "^^ms + + • • • + O J || £2 (£>.^- 1 ) 


— Il'^snap 


^ms 

l|2 


z;2(D;k-i) + l|ai</>i H + «J<('j|l£2(n;K-i) 


+ 2 / K (r^snap ~ l^ms) ' + ’ ’ ' + Oij4>j) 


Id 


— L^sna.n 


b'snap 
- 2 


|£2(D;k- 1) + ll«l</’l "I + ^j4‘j\\c'^{D-k,-^) 

' H- 

Jd 

Recall that Ui, ..., Qj are non-overlapping and each (j)j is supported on Qj. Therefore, 

||q!i0i -I- • • • -I- cr — 2k • (ai(/>i -I- • • • -|- aj(f)j) 

J D 

= -(l!b-oJiy +--- + \\RnAl. ). 


Hence, we have 


U^snaD 


,m+l II 2 


-1^ ^ 11‘^snaD '^^1 


\C^{D-,k-'^) ^ ll^snap 


= 1 - 


,m ||2 


.7 

msllZ;2(D;K-i) “ F/ 11-^^. 

1 = 1 




Z)Li II 


'i"y* 


lUsnap I’msil r2i'n-,^-i1 


'^snap 




Z)Li ll-^fij 


= 1 - 


\* snap J 


msll£2(£|.„-l) 


lllsnap llmsllc^^^j.^-ip 


by Lemma Hence the proof is complete. □ 

In Theorem we can see that the convergence rate of the online adaptive method can be small if the 
term is large. This term is determined at the beginning of the method when the initial number 

(i) 

of basis functions U is fixed. Therefore, we should choose h so that A) is significantly large. We will 
demonstrate the effect of choosing different initial number of basis functions in the next section. 


4 Numerical Results 

In this section, we will present some examples using both the offline and the online adaptive methods. To 
test the efficiency of the methods, we will compare the error of the solution from the adaptive methods with 
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the error of the solution obtained by uniform enrichment of the space, i.e. we increase the number of basis 
functions in all the coarse grid neighborhoods uniformly. We will also compare the error of the solution using 
the online adaptive method with different initial number of basis functions on each coarse grid neighborhood 
and different contrast in the permeability field. From the comparison in the online adaptive method, we 
can see that the convergence rate depends on the initial number of offline bases used. And if the initial 
basis is chosen in the appropriate way, the convergence of online adaptive method will be independent of the 
contrast. In the examples, we will use the following permeability fields with background value one. 



(a) Ki (b) K 2 


Figure 1: Permeability fields with high contrasts (denoted in red). 

Figure shows the permeability fields with background value one (shown in blue) and high contrasts 
(shown in red). In the numerical examples, we will vary the contrast so as to see the convergence rate of 
the adaptive methods with different contrast values. We will use the follow snapshot error to indicate the 
accuracy of the methods. 

^ ll^isnap Ums\\c^(D-,K~^) 

ll'^snap IIc2(D;k-1) 

4.1 Comparing the adaptive methods with uniform enrichment 

We compare the efficiency of the adaptive methods in this example. Consider equation 0 on the domain 
[0, 1]^ with homogeneous boundary condition, i.e. g = 0. We use coarse grid size 15 x 15 and fine grid 
size 40 X 40 on each coarse grid. We use the permeability field ki with contrast values le4 and le-4. The 
source function / is set to be 1 on top left coarse grid block, —1 on bottom right coarse grid block and zero 
elsewhere. We form the snapshot basis using spectral problem 1. 

We solve the equation in the following ways. 

Offline adaptive method: We use the offline adaptive method with initial number of bases per coarse 
grid neighborhood equal to 2, 0 = 0.2 and Jg = 0.5. 

Online adaptive method: We use the online adaptive method with initial basis functions obtained from 
spectral problem 1. The regions Hi,..., are 

(a) non-overlapping coarse grid neighborhoods; and 

(b) non-overlapping 2x2 coarse grid blocks. 

Uniform enrichment: we solve the equation again by increasing the number of bases in each coarse grid 
neighborhood uniformly from 2 to 40. 


15 





The results are compared by plotting the error e against the total number of bases used (Figure [^. 




(a) contrast = le4 (b) contrast = le-4 

Figure 2: Snapshot error of offline and online adaptive methods compared to the snapshot error of uniform 
enrichment for contrast le4 and le-4. 


DOF 

e (uniform) 

e (offline) 

2(840) 

0.0452 

0.0452 

8(3360) 

0.0115 

0.0017 

14(5880) 

0.0059 

2.75e-4 

20(8400) 

0.0039 

8.49e-5 

26(10920) 

0.0019 

2.98e-5 

32(13440) 

9.57e-4 

7.54e-6 

38(15960) 

1.60e-4 

7.84e-7 


DOF 

e (uniform) 

e (offline) 

2(840) 

0.0567 

0.0567 

8(3360) 

0.0123 

0.0012 

14(5880) 

0.0066 

2.71e-4 

20(8400) 

0.0040 

7.57e-5 

26(10920) 

0.0023 

2.11e-5 

32(13440) 

0.0011 

6.Ole-6 

38(15960) 

5.76e-4 

4.61e-7 


(a) contrast = le4 (b) contrast = le-4 

Table 1: Snapshot error of the offline adaptive method compared to the GMsFEM with the same number of 
basis functions. 


DOF 

e (a) 

e (b) 

2(840) 

0.0452 

0.0452 

3(1260) 

0.0056 

0.0045 

4(1680) 

0.0028 

1.36e-4 

5(2100) 

4.56e-4 

5.19e-6 

6(2520) 

9.57e-5 

4.94e-7 

7(2940) 

4.97e-6 

3.22e-9 

8(3360) 

1.03e-7 

1.14e-ll 

9(3780) 

4.84e-10 

6.80e-13 


DOF 

e (a) 

e (b) 

2(840) 

0.0567 

0.0567 

3(1260) 

0.0065 

0.0042 

4(1680) 

0.0033 

8.71e-5 

5(2100) 

6.20e-4 

4.08e-6 

6(2520) 

2.56e-4 

1.74e-7 

7(2940) 

9.21e-6 

2.81e-9 

8(3360) 

1.90e-7 

3.42e-ll 

9(3780) 

4.73e-10 

5.58e-12 


(a) contrast = le4 (b) contrast = le-4 

Table 2: Snapshot error of the online adaptive method using two different choices of Hi,..., llj. 

From Figure we can see that the snapshot error from the adaptive methods is always smaller than 
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the snapshot error obtained from uniformly increasing the number of basis functions. Moreover, the rate of 
convergence of the online adaptive method is faster than that of the offline adaptive method. Note that the 
online approaches require computations during the online stage of the simulations. From Table we can 
also see the difference in the convergence rate between the offline adaptive method and uniform enrichment. 
This shows that the error indicator in the offline method can successfully show the coarse grid neighborhoods 
with insufficient bases. 

For the online adaptive method with the two choices of regions, the one with 2x2 coarse grid blocks 
give a faster convergence rate, as observed from the Table This can be explained by the number of coarse 
grid neighborhoods a region contains. In (a), each region contains one coarse grid neighborhood while in 
(b), each region contains four. Therefore, the space Vh; for the calculation of the projection (j)j is larger in 
(b) than in (a) and captures more distant effects. Hence, the result is better in (b). 

4.2 Online adaptive method with different number of initial basis functions 

In this example, we focus on the online adaptive method and we want to see the effect of using different 
number of initial basis functions in the method. We consider permeability fields K 2 with the contrast le-4. 
For We divide the domain [0, 1]^ into 8x8 coarse grids and divide each coarse grid into 32 x 32 fine grids. 
The source function / is the same as the previous example. In each enrichment level, the regions Hi,..., Hj 
are chosen to be disjoint coarse grid neighborhoods. 

We solve the equation using 1, 2, 3 and 4 initial basis functions obtained from solving the two spectral 
problems. We plot the snapshot error e against the number of basis functions used in Figure|^and the value 
of e is shown in Tables [3] and ID 




(a) spectral problem 1 (b) spectral problem 2 

Figure 3: Snapshot error of online adaptive method with different number of initial basis functions. 

From Figure one can observe that if we use only 1 initial basis function, the rate of convergence is slow 
at the beginning. Similar behaviour can be seen if 2 initial basis functions are used, yet the convergence is 
faster. Using 3 initial basis functions seems to be the optimal choice in the sense that the snapshot error 
cannot be smaller when more initial basis functions are used. This can be explained by the value Amin. 
In the online adaptive method, this value depends on the initial basis functions obtained from the spectral 
problems. Theorem shows that the rate of convergence is bounded above by a value which decreases 
when Amin increases. When spectral problem I is used, the values of Amin are 0.0093, 0.0146, 2.5183 and 
5.0068 when 1, 2, 3 and 4 initial basis functions are used respectively. When spectral problem 2 is used, 
the values of Amin are 0.0016, 0.0247, 0.4939 and 0.7881 when 1, 2, 3 and 4 initial basis functions are used 
respectively. Therefore, the values of Amin corresponding to the first two basis functions are small. This 
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DOF 

e (1 basis) 

e (2 basis) 

e (3 basis) 

e (4 basis) 

1(112) 

0.2575 

/ 

/ 

/ 

2(224) 

0.1454 

0.0928 

/ 

/ 

3(336) 

0.1254 

0.0371 

0.0357 

/ 

4(448) 

0.0980 

0.0326 

0.0062 

0.0211 

5(560) 

0.0887 

0.0200 

0.0023 

0.0031 

6(672) 

0.0880 

0.0042 

1.68e-4 

7.29e-4 

7(784) 

0.0737 

1.73e-4 

1.68e-5 

6.02e-5 

8(896) 

0.0269 

8.65e-7 

1.87e-7 

2.56e-6 

9(1008) 

0.0019 

3.03e-9 

1.99e-9 

2.04e-8 


Table 3: Snapshot error of the online adaptive method using spectral problem 1 with 1 initial bases (Amin = 
0.0093), 2 initial bases (Amin = 0.0146), 3 initial bases (Amin = 2.5183) and 4 initial bases (Amin = 5.0668). 
Contrast is le-4. 


DOF 

e (1 basis) 

e (2 basis) 

e (3 basis) 

e (4 basis) 

1(112) 

0.5237 

/ 

/ 

/ 

2(224) 

0.1684 

0.2616 

/ 

/ 

3(336) 

0.1417 

0.0612 

0.1215 

/ 

4(448) 

0.1319 

0.0498 

0.0116 

0.0625 

5(560) 

0.1293 

0.0098 

7.71e-4 

0.0017 

6(672) 

0.1291 

2.48e-4 

1.46e-5 

3.49e-5 

7(784) 

0.0954 

4.39e-6 

1.84e-7 

4.10e-7 

8(896) 

0.0331 

2.51e-8 

5.31e-10 

1.08e-9 

9(1008) 

0.0034 

4.25e-ll 

1.83e-12 

1.17e-12 


Table 4: Snapshot error of the online adaptive method using spectral problem 2 with 1 initial bases (Amin = 
0.0016), 2 initial bases (Amin = 0.0247), 3 initial bases (Amin = 0.4939) and 4 initial bases (Amin = 0.7881). 
Contrast is le-4. 


suggests a criterion for choosing the initial number of basis function which is to include all basis functions 
with small eigenvalue from the spectral problem. 

Remark that the magnitude of Amin depends on the choice of the spectral problem. For spectral problem 
1, both the constant Cerr and Amin grow with the ratio H/h. For spectral problem 2, Cerr is independent of 
the mesh size and Amin is always bounded above by 1. 

4.3 Online adaptive method with different contrasts 

Next, we want to see the effect of varying the contrast to the online adaptive method. Similar to the previous 
example, we will start with different number of initial basis functions. Equation Q is solved in permeability 
field Ki with three different contrasts le-2, le-4 and le-6. The coarse grid size is 15 x 15 and the fine grid 
size is 40 x 40 on each coarse grid. The source function / is, again, 1 at top left corner and —1 at bottom 
right corner. For each number of initial basis functions, we plot the snapshot error of the method with the 
three different contrasts (Figures and [^. We also list the value of the snapshot error in Tables and 

From Figure]^ we can see that when spectral problem 1 is used, the change in the contrast has almost no 
effect on the snapshot error if we start with 2, 3 or 4 basis functions on each coarse grid neighborhood. This 
can also be confirmed by looking at Table However, if we use only 1 initial basis function, the contrast 
makes a huge difference. The convergence rate decreases as the contrast changes from le-2 to le-4 and then 
le-6. This result can also be explained by the value of Amin- The captions in Tablelist the values of Amin 
corresponding to different number of initial basis functions and different contrast value. These values are 
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(c) 3 initial basis functions. (d) 4 initial basis functions. 

Figure 4: Snapshot error of online adaptive method using spectral problem 1 with different number of initial 
basis functions and different contrast. 


consistent with Figure]^ since we can observe almost no changes in A„iin as the contrast varies if we use 2, 
3 or 4 initial basis functions, while decreases with the contrast for the case of 1 initial basis function. 

Similar behaviour can be observed when spectral problem 2 is used. We can observe from Table that 
if we start with one basis function per coarse grid neighborhood, the larger the contrast, the slower is the 
convergence rate. However, Figure also suggests that spectral problem 2 is less resistent to changes in the 
contrast than spectral problem 1. When 2, 3 or 4 initial basis functions are used, the snapshot error is nearly 
the same for contrast le-2 and le-4, yet the error increases as the contrast changes from le-4 to le-6. This 
jump is larger for 4 initial basis functions than 2 initial basis functions. Although spectral problem 2 is not 
as good as spectral problem 1 in the aspect of the resistance to changes in contrast, by comparing Tables 
and|^ we can see that the convergence rate is faster using spectral problem 2 than spectral problem 1. 

From the above observations, we know that the convergence rate of the online adaptive method can be 
independent of the contrast if we choose the initial basis functions well. Basis functions corresponding to 
small value of Amin are contrast-dependent, and we should include them in the initial basis. When Amin is 
large, we may expect that the online adaptive method will perform well with large contrast. 
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DOF 

le-2 

le-4 

le-6 

1(420) 

0.1772 

0.2039 

0.2130 

2(840) 

0.0560 

0.0753 

0.0839 

3(1260) 

0.0399 

0.0620 

0.0699 

4(1680) 

0.0151 

0.0597 

0.0676 

5(2100) 

0.0039 

0.0595 

0.0674 

6(2520) 

3.73e-4 

0.0567 

0.0674 

7(2940) 

1.47e-4 

0.0317 

0.0643 

8(3360) 

1.22e-7 

0.0010 

0.0466 

9(3780) 

2.43e-10 

8.41e-6 

0.0035 


(a) 1 initial basis function. Values of Amin are 
0.0776, 0.0027 and 5.12e-5 for contrasts le-2, 
le-4 and le-6, respectively. 


DOF 

le-2 

le-4 

le-6 

3(1260) 

0.0341 

0.0317 

0.0303 

4(1680) 

0.0041 

0.0037 

0.0037 

5(2100) 

0.0016 

0.0014 

0.0014 

6(2520) 

2.40e-4 

2.15e-4 

2.02e-4 

7(2940) 

3.26e-5 

2.74e-5 

2.19e-5 

8(3360) 

7.28e-7 

6.45e-7 

6.13e-7 

9(3780) 

1.15e-8 

3.86e-9 

6.83e-9 


(c) 3 initial basis functions. Values of Amin are 
3.6204, 3.5980 and 3.5978 for contrasts le-2, 
le-4 and le-6, respectively. 


DOF 

le-2 

le-4 

le-6 

2(840) 

0.0472 

0.0452 

0.0399 

3(1260) 

0.0060 

0.0056 

0.0054 

4(1680) 

0.0030 

0.0028 

0.0027 

5(2100) 

5.47e-4 

4.56e-4 

4.42e-4 

6(2520) 

1.41e-4 

9.57e-5 

9.04e-5 

7(2940) 

l.lle-5 

4.97e-6 

5.29e-6 

8(3360) 

1.20e-7 

1.03e-7 

l.lOe-7 

9(3780) 

2.67e-10 

4.84e-10 

6.16e-10 


(b) 2 initial basis functions. Values of Amin are 
1.9511, 1.9264 and 1.9262 for contrasts le-2, 
le-4 and le-6, respectively. 


DOF 

le-2 

le-4 

le-6 

4(1680) 

0.0257 

0.0240 

0.0239 

5(2100) 

0.0024 

0.0024 

0.0024 

6(2520) 

6.92e-4 

6.98e-4 

6.99e-4 

7(2940) 

9.78e-5 

1.Ole-4 

9.92e-5 

8(3360) 

7.32e-6 

6.63e-6 

6.62e-6 

9(3780) 

1.26e-7 

9.01e-8 

8.12e-8 


(d) 4 initial basis functions. Values of Amin are 
5.2765, 5.2396 and 5.2656 for contrasts le-2, 
le-4 and le-6, respectively. 


Table 5: Snapshot error of online adaptive method using spectral problem 1 with different number of initial 
basis functions and different contrast. 


5 Conclusion 

In this paper, we present two adaptive enrichment algorithms for the mixed GMsFEM. The first one is an 
offline adaptive method which adds basis computed from the offline stage. This method is based on a local 
error indicator which is the norm of a residual operator restricted on a local space. Offline basis functions 
are added to those coarse grid neighborhoods with large errors. The other algorithm is an online adaptive 
method. Online basis functions are constructed by solving local problems based on a residual operator. We 
select non-overlapping regions in the domain and solve for an online basis function on each of the regions. 

We show theoretically the convergence of the two methods. For the offline method, the rate of convergence 
depends on two parameters, which control the number of coarse grid neighborhoods to be selected to add 
basis functions and the number of basis functions to be added. The convergence rate of the online method 
depends on the eigenvalues corresponding to those basis functions that are not included in the initial basis. 
The larger the above mentioned eigenvalues, the faster is the convergence rate. The numerical results are 
consistent with these hndings. It is also shown numerically that if the basis functions corresponding to the 
eigenvalues that are contrast-dependent are included in the initial basis, the online adaptive method will 
be resistant to the change in the contrast value. Those eigenvalues are the smallest ones and therefore the 
corresponding basis functions should be included in the initial basis to speed up the convergence. 
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(c) 3 initial basis functions. (d) 4 initial basis functions. 

Figure 5: Snapshot error of online adaptive method using spectral problem 2 with different number of initial 
basis functions and different contrast. 
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DOF 

le-2 

le-4 

le-6 

1(420) 

0.5506 

0.5716 

0.5810 

2(840) 

0.0799 

0.0968 

0.1054 

3(1260) 

0.0477 

0.0719 

0.0803 

4(1680) 

0.0021 

0.0692 

0.0778 

5(2100) 

0.0045 

0.0690 

0.0776 

6(2520) 

9.73e-4 

0.0649 

0.0776 

7(2940) 

2.23e-5 

0.0463 

0.0735 

8(3360) 

2.20e-7 

0.0298 

0.0401 

9(3780) 

3.42e-10 

4.1024e-4 

0.0184 


(a) 1 initial basis function. Values of Amin are 
0.0510, 9.41e-4 and 4.33e-5 for le-2, le-4 and 
le-6 respectively. 


DOF 

le-2 

le-4 

le-6 

3(1260) 

0.1202 

0.1418 

0.1319 

4(1680) 

0.0128 

0.0148 

0.0357 

5(2100) 

0.0014 

0.0020 

0.0043 

6(2520) 

6.66e-5 

7.69e-5 

5.68e-4 

7(2940) 

1.02e-6 

1.73e-6 

1.14e-5 

8(3360) 

5.19e-9 

9.20e-9 

9.76e-8 

9(3780) 

9.61e-12 

2.29e-ll 

1.13e-9 


(c) 3 initial basis functions. Values of Amin are 
0.6460, 0.4882 and 0.5763 for le-2, le-4 and 
le-6 respectively. 


DOF 

le-2 

le-4 

le-6 

2(840) 

0.2909 

0.3057 

0.3164 

3(1260) 

0.0270 

0.0283 

0.0477 

4(1680) 

0.0050 

0.0052 

0.0116 

5(2100) 

6.66e-4 

9.15e-4 

0.0017 

6(2520) 

3.44e-5 

4.67e-5 

8.54e-5 

7(2940) 

8.12e-7 

1.02e-6 

2.53e-6 

8(3360) 

4.13e-9 

5.20e-9 

2.06e-8 

9(3780) 

1.76e-ll 

2.11e-ll 

8.78e-ll 


(b) 2 initial basis functions. Values of Amin are 
0.3176, 0.3823 and 0.3790 for le-2, le-4 and 
le-6 respectively. 


DOF 

le-2 

le-4 

le-6 

4(1680) 

0.0824 

0.0930 

0.1128 

5(2100) 

0.0046 

0.0044 

0.0315 

6(2520) 

1.84e-4 

1.99e-4 

0.0035 

7(2940) 

2.17e-6 

3.15e-6 

4.84e-5 

8(3360) 

1.05e-8 

1.23e-8 

3.25e-7 

9(3780) 

2.02e-ll 

3.65e-ll 

1.02e-9 


(d) 4 initial basis functions. Values of Amin are 
0.7697, 0.7650 and 0.7307 for le-2, le-4 and 
le-6 respectively. 


Table 6: Snapshot error of online adaptive method using spectral problem 2 with different number of initial 
basis functions and different contrast. 
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